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Abstract. The cosmic web contains a large fraction of the total gas mass in the universe but is difficult to detect at most 
wavelengths. Synchrotron emission from shock-accelerated electrons may offer the chance of imaging the cosmic web at radio 
wavelengths. In this work we use 3D cosmological ENZO-MHD simulations (combined with a post-processing renormalisation 
of the magnetic field to bracket for missing physical ingredients and resolution effects) to produce models of the radio emission 
from the cosmic web. In post-processing we study the capabilities of 13 large radio surveys to detect this emission. We find that 
surveys by LOFAR, SKA 1-LOW and MWA have a chance of detecting the cosmic web, provided that the magnetisation level 
of the tenuous medium in filaments is of the order of ~ 1% of the thermal gas energy. 
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1. Introduction 


In cosmological structure formation, tiny matter perturba¬ 
tions grow through supersonic accretion of smooth and cold 
(T < 10 4 K) matter. In the process, matter gets shock heated, 
thereby converting the infall kinetic energy into thermal 
energy ( Sunvaev & Zeldovichl Il972h . Cosmological simula¬ 


tions have been used to measure the distribution of Mach 
numbers in cosmic shocks using a variety of methods. It has 
been found that weak (M ^2 — 4) merger shocks inside 
galaxy cluster are responsible for most of the conversion of 
the kinetic energy into gas heating, while stronger external 


ripheral regions of large-scale structures (e.g. Rvu et al. 

2003; 

Pfrom 

mer et al. 2006; Skillman et al. 2008; 

Vazza et al. 

2009, 

201 i 

Arava-Melo et al. 2012c Schaal & Springel 2014 

). All 


mass in the Warm-Hot-Intergalactic-Medium (WHIM) with 
temperatures in the range ~ 10 5 — 10 7 K has been processed 
by accretion shocks. Up to « 90 percent of the total baryon 
and dark matte r mass in the Universe is expected to r eside in 
this gas phase dCen & Ostriked[l999l:iDave et al.ll200ll) . which 
is extremely difficult to be observed at most wavelengths. 
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Indeed, only a few tentative detections of the WHIM associ- 
ated with cosmic filame n ts have been c laime d in X-ray (e.g. 


Finog uenoy et al.1 120031: IWerner et all 120081: iNicastro et al 


2010L12013) or at microwaves using the Sunyaev Zeldovich 


effect (Planck Collaboration et al. 2013). A few possible 


detections of s hocks emitting in radio have b een also reported 
( Bagchietal. 2002 : Kronberg et al.1 120071: iGiovannini et al.1 
2010 : iFarns worth et al.11201 3b . It has been pointed out that 


the accretion shocks around cosmological filaments may 
be detected with radio telescopes provided these shocks 
are efficie nt enough in acc e lerating electrons to relativistic 


energies (Keshetetal. 2004; Brown 2011; Arava-Melo et al 


2012 : IVazza et al.1 2015bt IGiovannini et al. lio 1 5h . A few 
authors have studied the emission from electrons accelerated 
at th e intersections of galaxy clusters and cosmic filaments 


(e.g. Keshet et al.112004 iBattaglia et al.1120091: ISkillman et al 


2011), yet the detections were mostly related to intermittent 


mergers of g alaxy clusters, rather than to stationary accretion 
shocks. iBagchi et al.l d2002h used the numerical method of 
iMiniati et al l ( 2001 ) to simulate a scenario in which a filament 
shock causes the observed radio emission in the cluster ZwCl 
2341.1+0000. In order to explain the observed radio emission, 
~ 2 percent of the shock ram energy has to go into relativistic 
electrons, and a magnetic energy of ^ 1.5 fiG is required, 
probably indicating that the emission is related to more stan¬ 
dard merger shocks dvan Weeren et al.ll2009l: IGiovannini et al.1 
201 Oh . 
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Table 1 . List of the simulations run for this project. First col¬ 
umn: box length the simulated volume; second column: num¬ 
ber of grid cells in the initial conditions; third column: spa¬ 
tial resolution; fourth column: DM mass resolution; fifth col¬ 
umn: note on the adopted physics (MHD or Cosmic Rays (see 
Appendix). 


T'box 

N grid 

Ax 

muM 

note 

[Mpc/h] 


[kpc/h\ 

[Mo/h] 


50 

2400 3 

20.8 

9.8 • 10 b 

MHD 

100 

1200 3 

83.3 

6.2 • 10 7 

MHD 

200 

1200 3 

166.6 

4.9 • 10 8 

MHD 

300 

2048 3 

146.5 

3.4 • 10 8 

CR 


More recently, lArava-Melo et al.1 ( 20121) studied the 
emission from primary electrons accelerated at cosmological 
shock waves and concluded that filaments must host fewer 
radio objects than galaxy clusters, with an average of a few 
radio objects brighter than 10 32 erg/(s • Hz) at 1.4 GHz within 
(500 Mpc/h) 3 , and that most of this emission should come 
from the intersections with surrounding galaxy clusters. They 
predicted a flux density of radio emission from filaments of 
the order of 0.12 /iJy at redshift z ~ 0.15 and at a frequency 
of 150 MHz. 


We recently found that the detection of the biggest fila¬ 
ments of the cosmic web might be within reach of the Square 
Kilometer Array if the magnetic field there reaches ^ 10 — 
100 nG dVazza et al.1 l2015bb . While these first results were 
based on the simple assumption of a fixed magnetic field for 
a representative filament of our sample, in this work we use 
MHD simulations with ENZO ( Brvan et al. 2014 ) and study a 
much larger volume. 


2. Methods 

2.1. Numerical simulations with enzo 


ENZO is a highly parallel code for cosmological simulations, 
which uses a particle-mesh N-body method (PM) to follow the 
dynamics of the DM and offers a variety of MHD-hydro solvers 
to compute the evolution of cosmic gas. 

Here we mostly anal yze simulations obtained with ENZO- 
MHD described in IVazza et aD ( 2014a ), which belong to a 
large project run on Piz-Daint (Lugano, Switzerland) for a 
CHRONOS 0 proposal. This version of ENZO isj^ased on the 
Dedner formulation of MHD equations dDedner et al.112002 ). 
which uses hyperbolic divergence cleaning to keep the V • 
B as small as possible during the computation. The MHD 
solver adopted here uses a piecewise-linear method recon¬ 
struction (PLM), where fluxes at cell interfaces are calcu¬ 
lated using the local Lax-Friedrichs Riemann solver (LLF 
Kurganov & TadmodhOQOl) . Time integration is performed us¬ 
ing a total variation diminishing (TVD) second order Runge- 
Kutta (RK) scheme dShu & Oshd 1988). The PLM+MHD 
solver used in the simulations presented here (as well as the 


http://www.cscs.ch/user_lab/allocation_schemes/index.html 


non-MHD version of the PPM hydro solver), has been ported to 
Nvidia's CUDA framework, allo wing ENZO to take advantage 
of mo dern graphics hardware dWang et al.1 l2010t iBrvan et all 


20141) . In simulations using a single-level mesh, the GPUs pro¬ 
vide a factor of ~ 4 in speed-up over CPU-only runs. Our run 
with the highest spatial resolution is a (50 Mpc) 3 volume sim¬ 
ulated using 2400 3 cells and 2400 3 DM particles (resolution 
20.8 kpc), which, as far as we know, is t he la rgest MHD cosmo¬ 
logical simulation to date dVazza et al.ll2014ah . This run used 
^4.5 million core hours running on 512 nodes (2048 cores in 
total) on Piz Daint in Lugano. Fig.Q]shows the projected (mass- 
weighted) temperature and magnetic field strength at z = 0 for 
our best-resolved run (2400 3 ). The intricate network shows the 
outer shell of accretion shocks where kinetic energy is ther- 
malised. In the inner regions, in the centres of halos and along 
the spine of filaments, the magnetic field gets amplified. A list 
of the simulations analysed for this work is given in Tab. 1. To 
extend our results to a larger volume and include the role of sec¬ 
ondary electrons injected by cosmic ray protons, we analyse a 
set of simulations that include CRs (details in the appendix and 


in 


Vazza et aDll2Q14bb . 


2.2. Cosmological Model 

We assume a WMAP 7-year cosmology dKomatsu et al.ll20lH) 
with Do = 1.0, Db = 0.0455, Ddm = 0.2265, Da = 0.728, 
Hubble parameter h = 0.702, a normalisation for the primor¬ 
dial density power spectrum as = 0.81 and a spectral index of 
n s = 0.961 for the primordial spectrum of initial matter fluc¬ 
tuations. All runs start at Z[ n > 50, the exact value of which 
varies with resolution. 


2.3. Non-thermal processes 

In order to model the synchrotron emission we have to know 
the energy spectrum of relativistic electrons and the magnetic 
field at each point in space and time. Both are largely unknown, 
specially at the low-density regime of cosmic filaments. In our 
analysis, we try to bracket the possible range by varying model 
parameters. 


2.3.1. Magnetic fields 

The degree of magnetisation in the low-density cosmic web 
is largely uncertain. The predictions from cosmological sim¬ 
ulations lie in the range of 10 4 /iG — 0.1/iG for fila¬ 
ments, depending on numerical resolution of the cosmolog- 


nario (e.g. Dolag et al. 1999; Br 

ti 

ggenetal. 2005; Sig 

et al. 

2003; Rvu et a 

2008; Xu et al. 

2009; Donnertetal. 

o 

o 

<N 

Marinacci et al. 

2015). Since the growth of primordial mag- 


netic fields in filaments is do minated by compressive mo 
tions and small-scale shocks dRvu et al. 120081: IVazza et al 


2014al) . the dynamical memory of the system might persist over 
long cosmological times and the observed magnetisation level 
closely connects to the seeding process(es). This is different 
from galaxy clusters where most of the magnetic energy is ex- 
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Fig. 1. 3D rendering of the projected temperature (purple colours for T < 10 5 K and red for T > 10 5 K) and magnetic field 
intensity (yellow for B > lOnG and orange for B < lOnG) for our simulated (50Mpc) 3 volume at z = 0, simulated with 2400 3 
cells and DM particles. 


tracted from the kinetic energy, thereby quickly erasing previ¬ 
ous dynamical information. In our latest MHD runs we investi¬ 
gated the growth of magnetic fields startin g fro m B o = 10~ 1Q 
G (comoving), seeded at high redshift (e.g. Widrowet alJ201ll 
iRvu et al.1 1201 ll and references therein). However, the mag 
netic field measured even in our most massive clusters hardly 
reaches ^0.01 — 0.05 11 G, while in filaments the magnetic field 
is never larger than a few nG. In clusters, the above values are ^ 
10 — 100 times smaller than the typical magneticfield inferred 
using the Faraday Rotatio n effect (e.g. iMurgia et al.1 12004 : 


Bonafede et al .201(1 20131) . Therefore, something appear to be 


missing in our numerical modeling of magnetic fields. On the 
one hand, only within galaxy clusters we might have a large 
enough dynamical range t o simulate a small-sca le dynamo(e.g. 
Schekochihin et al. 20 q 4 I Schleicher et akll2013l) . In filaments, 


however, the dominance of compressive forcing makes dy¬ 
namo amplification inefficient. Even in high-resolution adap¬ 
tive mesh refinement runs we did not measure any strong field 
growth beyond that caused by mere compression and mag¬ 
netic flux conservation. On the other hand, additional seed¬ 
ing from astrophysical sources (star forming regions, galax¬ 
ies and galaxy clusters) can magnetise large-scales, even if 


1999; Donnert et al. 2009; Xu et al. 2009; Schober et al. 

2013; 

Beck et al. 2012)7Fast-growing plasma in 

stabilities, such as 
/3 plasmas (where 
; magnetic energy) 
p to values we ob- 

firehose and mirror instabilities, in large 
/3 is the ratio between the thermal and the 
might be able to amplify magnetic fields u 

serve in clusters (Schekochihin et al. 2005 

; Kunz et al. 2011; 


tional seeding for the dynamo. Indications of substantial mag¬ 
netisation along the filamentary region in the SW sector of 
the Coma cluster have been reported in lBonafede et al.1 ( 2013 ). 
Finally, it has also been suggested that cosmic-ray particles 
accelerated at stron g shocks can amplify upstream magnetic 
fields substantially dDrurv & Downesll2012 : Bruggen|[2013l) . 


In summary, there is no firm conclusion about the strength 
of magnetic fields outside galaxy clusters, neither from obser¬ 
vations nor from theory. Future radio observations will be cru¬ 
cial to assess this, and our goal here is to provide quantitative 
predictions for different magnetic field models, starting from 
extreme cases. Hence, in post-processing we renormalize the 
magnetic field using two separate recipes: a high-amplification 
model (HA) and a low-amplification model (FA). In both cases, 
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Fig. 3. The £(M, T) function we use to compute the acceler¬ 
ation of CR-electrons at shocks, based on iHoeft & Bruggenl 
d2007t) . 


the normalisation depends on the local gas overdensity, n/n cv , 
where n is the gas density and n cr?g = 3Hq /(87rG)£2b/^M 
is the critical density rescaled for the gas fraction. In the HA 
case, we model the efficient magnetisation of all large-scale 
structures and increase the magnetic field for n/n cr > 2. In 
the LA case, we limit the renormalisation to densities larger 
than the gas density at the virial radius of resolved halos in 
the simulation, n/n CY > 50. In both cases, the strength of 
all components of the magnetic field in each cell is renor¬ 
malized to (3 = 100 , where (3 is the ratio between the ther¬ 
mal energy and the magnetic energy. On the scale of clusters 
and filaments, the HA model typically yields a ^ 10 times 
larger magnetic field compared to what is predicted by our 
MHD runs, up to several ~ /iG. The projected maps of mag¬ 
netic fields and the (F?, n/n cr ) phase diagrams for these two 
models are given in Fig. [2] In the density range of filaments 
(n ~ 5 — 30n cr ) the prediction of the two models differ sig¬ 
nificantly, i.e. B rsj 0.001 fiG in the LA model (equivalent to 
our original MHD run for n < 50n cr ) and B ~ 0.05 /iG in the 
HA model. For B < 3.2 /iG and flat electron spectra the syn¬ 
chrotron emission scales as P v oc B 2 (EqUbelow) and, hence, 
this difference can have a large impact on the detectability of 
filaments in radio. Cells in the regime n ~ 10 — 100 n cr are 
more relevant for giant filaments and the outskirts of galaxy 
clusters, and the magnetic field varies from 0.001 — 0.01 /iG 
(MHD model) to 0.1 — 1 /iG (HA and LA). As we will show, 
the detectability of filaments in radio will crucially depend on 
what is the level of magnetisation at their typical overdensity, 
and therefore radio observation may have the potential of prob¬ 
ing the mechanisms for magnetic field amplification in these 
regions. 


2.3.2. The acceleration of relativistic particles at 
shocks 

A plethora of mechanisms can contribute to the accelera¬ 
tion of electrons, e.g., direct injection from galactic winds 
and AGN, acceleration by cosmological shock waves, contin¬ 
uous injection by hadronic collisions and turbulent reaccel¬ 
eration. Simulations suggest that a vast majority of relativis¬ 
tic electrons in the WH IM should be produced by cosmolog¬ 
ical shock waves (e.g. iRvu et alJl2003c IPfrommer et alJl 200 (x 


Arava-Melo et al .1 120 1 2h . even though el ectron accelerati on in 

high-/3 plasmas is poorly understood (e.g. lGuo et al.l2014 ). The 
theory of Diffusive Shock Acceleration (DSA) offers a theoreti¬ 
cal framework to estimate the fraction of the kinetic energy that 
shocks dissipate into the acceleration of cosmic rays, and links 
the slope of the produced momentum distribution to the shock 
Mach number (e.g. lKang & Ryull2013l and references therein). 
However, the microphysical details are largely unknown and 
the efficiency of this process, specially at the weak shocks 
(M < 5) in clusters, are n ot understood (e.g. I Vazza & Bruggenl 
201 4 IVazza et al.ll2015al and discussion therein). In order to 
predict the synchrotron emission from acceler ated electrons, 
we use the formalism bv lHoeft & Bruggenl d2007l) that assumes 
an exponential cutoff in the energy distribution of electrons, de¬ 
termined by the balance of the acceleration rate from DSA and 
of the (synchrotron and Inverse Compton) cooling rate. In the 
downstream region, DSA is assumed to generate supra-thermal 
electrons that follow a power-law in energy. The total emission 
in the downstream region is computed by summing up the con¬ 
tributions from all electrons accelerated at the shock. 

The monochromatic radio power at frequency v, P u , de¬ 
pends on the shock surface area, S, the downstream electron 
density, n<\, the electron acceleration efficiency, £(M, xjl , the 
downstream electron temperature, Xd, the spectral index of the 
radio emission, s = <5/2, the magnetic field to the CMB energy 
density, Bqmb and the magnetic field in each cell, B : 


P v oc S ■ n d • £(M, T) ■ • T c 


3/2 


fil+5/2 


(^CMB + ^ 2 ) 


( 1 ) 


Relic relativistic electrons resulting from previous shocks are 
expected to accumulate around 7 ^ 10 2 dBrunetti et al.ll2QQlh 
and can be re-accelerated by subsequent shocks at much higher 
efficiency. The net effect of freshly injected and re -acce lerated 
particles can be modelled using the formalism bv lKang & Rvul 
( 2013 ). where the global effect of re-acceleration can be incor¬ 
porated into the acceleration efficiency of CRs, to correct the 
£(M, T) function of the HB07 model. It has been estimated in 
Pinzke et al.1 (120130 that for weak enough shocks in the ICM 
(M < 3) re-accelerated electrons dominate over directly in¬ 
jected electrons. For the sake of simplicity, we assume here that 
everywhere in the simulated volume the overall acceleration ef¬ 
ficiency of CRs for M <3 is affected by the re-acceleration of 
previously injected electrons. In this case the acceleration effi¬ 
ciency is ^ 10 -5 — 10 -6 for M = 3 (depending on the temper- 


2 We note that, in agreement with most of the literature on the DSA 
of electrons, we include in our definition of the electron acceleration 
efficiency the dependence on M in EqQ] i.e. £(M, T) = £ e ,o4i(M), 
where £ e ,o = 10 -3 and 4/M, T) is given in Eq. (31) of HB07. 
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LA model 



~z 


Iog10(n/n_cr) 


Fig. 2. Top: Mass-weighted magnetic field maps for the HA and LA model, for the same volume of Fig. [I] The colour-bar is in 
units of log 10 ([/xG]). Bottom: phase diagrams on the (L>, n/n cr ) plane for the same models. We additionally draw isocontours 
tenfold in the number of cells contributing to the pixels in the phase diagrams, starting from 100 cells. 


ature) and ~ 4 — 8 • 10 -4 for M = 10. Fig. Ogives the shape 
of the £(M, T ) function for different temperatures, and shows 
the role of shock re-acceleration for M < 10 shocks. Given 
the very steep dependence on Mach number of the acceleration 
efficiency at weak shocks, the presence of re-accelerated parti¬ 
cles can boost the net acceleration by one order of magnitude. 


Shock waves are identified in post-processing using the 
same procedure as in i Vazza et al.l (2009), which is based on the 
analysis of velocity jumps of cells preliminary tagged based on 
the local 3D divergence. In the case of multiple shocked cells 
along the shock propagation axis, we perform a further clean¬ 
ing in order to select only the cell yielding the largest Mach 
number. To mask out some small differences in the reioniza¬ 
tion background, in post-processing we imposed a temperature 
floor of T re = 5 • 10 3 K everywhere ( 1 Vazza et al.112009k . The up¬ 
stream/downstream shock parameters are then used in Eq.Q]to 


compute the synchrotron emission as a function of frequency 
and build full-sky emission maps of our simulated volumes. 
The procedure is run in parallel for both magnetic field models 
of Sec. 12.3.11 to produce a full catalog of mock radio observa¬ 
tions. 


Fig.[4]shows the emission-weighted projected map of mean 
radio spectral index along the line-of-sight, for our highest res¬ 
olution box and assuming the HA model for the magnetic field. 
The accretion shocks have flat radio spectra, s ~ 1 — 1.5, while 
weaker shocks within clusters or filaments have steeper spec¬ 
tra, _s_rZ_ 1J) — 3, as in the case of typical radio relics (e.g. 
Hoeft & Bruggenll2007l) . 
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Table 2. Fiducial parameters for our simulated radio observations. The quoted frequencies are reference values for each instru¬ 
ment. The last column gives the full resolution sensitivity we adopted in surface brightness, which is 3a of the assume thermal 
noise per beam (in the case of SKA1-LOW, this is the confusion noise). See Sec. l2.3.3l for details. 


array 

configuration/strategy 

frequency 

[MHz] 

resolution 

[arcsec] 

min. baseline 
[m] 

sensitivity 

[mJy/beam] 

detection threshold 

[/xJy/arcsec 2 ] 

VLA 

VLSSr 

74 

80 

35 

100.0 

42.365 

VLA 

NVSS 

1400 

45 

35 

0.45 

0.588 

Westerboork 

WENSS 

330 

54 

36 

3.6 

3.268 

GMRT 

TGSS 

150 

20 

100 

5.0 

33.098 

Molonglo 

SUMSS 

840 

43 

15 

1.0 

1.307 

LOFAR-HBA 

Large Survey 1 

120 

25 

40 

0.25 

1.059 

LOFAR-LBA 

Large Survey 

40 

25 

40 

2.0 

8.473 

LOFAR-HBA 

Large Survey 2 

120 

5 

40 

0.1 

10.591 

MWA 

Broadband Survey 

150 

120 

7.7 

10 

1.838 

SKA1-LOW 

Cont. Survey 

120 

10 

45 

0.02 

0.17 

SKA 1-MID 

Band2 Wide Survey 

1000 

0.5 

15 

0.001 

10.591 

SKA 1-MID 

Band2 Deep Survey 

1000 

0.5 

15 

0.0002 

2.118 

ASKAP 

EMU 

1400 

10 

12 

0.01 

0.264 



-4 - 3.7 - 3.4 - 3.2 - 2.9 - 2.6 - 2.3 -2 - 1.8 - 1.5 - 1.2 


Fig. 4. Emission-weighted projected mean radio spectral index 
for a subregion with side 27 Mpc (and depth 50 Mpc) through 
our highest resolution volume, assuming the HA model for the 
magnetic field. The additional contours show the projected gas 
pressure for the same volume. 


2.3.3. Mock radio observations 


In [Vazza et al. (2015b) we explored the possibilities of future 
surveys and deep pointings with the SKA to detect the radio 
emission from filaments in targeted observations. In general, 
predicting the emission from a specific filaments is not trivial 
and the spatial distribution of gas in filaments in the surround¬ 
ings of galaxy clusters is not known either. In this work, we 


instead explore the potential of current and future radio sur¬ 
veys to detect the diffuse radio emission from the cosmic web 
in a statistical way and without a priori knowledge of the best 
fields to target. Table 2 gives the full set of specification used 
to produce a mock observation of a simulated field, according 
to a procedure similar to IVazza et al. ( 2015b ): 


1. We compute the radio emission at different redshifts pro¬ 
duced during the simulation, in their comoving frame, and 
compute the total emission map by summing up the contri¬ 
bution of all cells along the line of sight. 

2. We convert the emitted power in the physical frame, as¬ 
suming an observer at different redshift, i.e. we include the 
effect of the luminosity distance of the emitting frame and 
the effect of cosmological dimming. 

3. We transform the images in the Fourier space using FFT. 
We remove from the Fourier plane the frequencies smaller 
than the minimum antenna baseline of each specific ra¬ 
dio configuration, mimicking the loss of signal from scales 
larger than the minimum instrumental baseline. This is par¬ 
ticularly relevant for the large-scale diffuse emission we are 
investigating here. 

4. After the baseline removal the images are converted back 
into real space. Then the emission is convolved for the 
beam using a 2D Gaussian filter. 

5. Only the pixels whose emission is I (x,y,z) > 3cr are con¬ 
sidered as detected by the specific observation. Depending 
on the array configuration, a is either the confusion limit 
(e.g. the contribution of unresolved background of contin¬ 
uum radio sources below the beam size) or the thermal 
noise of the instrument. The equivalent noise per beam is 
computed as 7 t hr ~ cr/(1.1333W / | eam ), where Wseam is 
the beam size. 


In all cases, we implicitly a ssume a perfect remova l of 


the Milky Way foreground (e.g. Bo naldi & Brownll2015L for 


the case of SKA observations), of all resolved point-like radio 
sources as well as an ideal calibration and deconvolution of the 
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relic emission in different volumes 



ProdJW/fg] 



10 21 10 22 10 23 1 0 24 10 25 10 26 10 27 
P r0 dio[W/Hz] 


Fig. 5. Radio distribution function for the clusters in our three 
simulated volumes at 1.4 GHz, using the LA model. No obser¬ 
vational effects are considered in this case. 


Fig. 6. Radio distribution function for the clusters in our 
(100 Mpc) 3 volume, normalized to Gpc 3 , for all investigated 
emission models. The additional gray lines show the range of 
the best fit to the distribution of radio relics based on the NVSS, 
computed by Nuza et al. (2012). 


radio data. We post-process our simulated boxes under many 
different observing configurations, exploring the capabilities 
of existing and future surveys in radio. The past and current 
surveys are: a) VLSS Redux survey at 74 MHz with the VLA 
(jLane et al.l2012h: b) NVSS survey at 1400 MHz with the VLA 


); c) the WENSS survey at 330 MHz with 
Westerbork ( Rengelink et all 19971): d) the TGSS survey at 150 
MHz with GMRT ( Bagchi et al. 2011 ): e) the SUMSS survey 
at 840 MHz with Molonglo ( Bock et al.lll999h . For the future 
surveys it is more difficult to quote specific numbers. Thus, in 
some cases we explored a few possible configurations bracket¬ 
ing their specifications. In detail, we studied: a) Large Area 
surveys with LOFAR HBA (120 MHz) and LBA (40 MHz) 
dRottgering et alJl^oTTh . The beam of LOFAR is formed dig¬ 
itally and the final resolution/sensitivity that each survey will 
be able to reach is not yet jixed. For this reason, we studies 
two cases, based on lCassano et al.l d20l4 : b) SKA1-LOW (120 
MHz) and SKA1-MID (1000 MHz) surveys. Also in this lat- 


ter case, we explored two 1 

possible SKA 1-MID configurations, 

as in Prandoni & Seymour 

(2014); c) Broad Band Survey with 

MWA at 150 MHz, as in 1 

ringav et al. (2013b: d) EMU survey 

with ASKAP at 1400 MHz (Norris et al.l 2013). 


All used parameters are given in Tab. 12.3.21 The last column 
gives the detection level in surface brightness for each survey. 

In the case of SKA1-LOW, this is computed considering confu¬ 
sion noise. We also list the field-of-view of each instrument but 
this never entered our calculations since we consider full-sky 
surveys or very large surveys obtained by adding many point¬ 
ings. 

A public repository of radio maps for the 

full volumes studied in this work is available at 
http:/Ayww.hs.uni-hamburg.de/DE/Ins/Per/Vazza/projects/PublicJfetb.(tanl>ines the exquisite uniform resolution of our 2400 3 

run with the large volume of the (200Mpc) 3 is presently un- 


3. Results 

3.1. Comparison with radio-relic emission in the 
NVSS 

First, we check our emission model against the available con¬ 
straints fr om the NVSS, us i ng rel ics in galaxy clusters. The 
model bv iHoeft & Bruggenl ( 2007 ) has been derived to model 
radio relics and in our case it fits rather well the distribution 
of relic emission in clusters. We compute the radio emission 
within i7 2 oo (~ 0.7 77 v i r ) of all identified clusters in the three 
simulated volumes at 2 = 0.05 and compare the observed 
statistics of radio relic emission for the same reference mag¬ 
netic field model (here the LA model), first without including 
observational effects in order to better focus on the effects of 
resolution and volume in our runs. The result is shown in Fig 13 
where we renormalise the measured statistics to the same vol¬ 
ume of (1 Gpc) 3 . Two clear trends stand out. First, the distri¬ 
bution functions of the different volumes have the same shape, 
with a factor of ~ 10 different normalisation in the number 
of objects when the high resolution and low resolution are 
compared. This follows from the fact that the internal dynam¬ 
ics (leading to merger shocks) of_ clusters is better resolved 
at higher resolution ( Vazza et al .1120 Ilk , producing more radio 
emitting objects. Second, the increase in the simulated volume 
enables the formation of more massive clusters, enabling more 
powerful major mergers. For this reason, our largest volume 
has a ^ 10 times more powerful tail of high radio power. 
The above trends are in line with what already known in the 


literature and based on shock statistics ( e.g. iRvu et al.l 12003 


Pfrommer et al.l 12006 ^ Vazza et al.l 2009112011 ) and on simu¬ 
lated radio relics dSkillman et al.l l^oTTh . Having a simulation 
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Fig. 7. Radio emission from a subregion with side 27 Mpc (and 
depth 50 Mpc) through our highest resolution volume for the 
HA and LA magnetic fields models (units of log 10 [Jy/arcsec 2 ] 
assuming a luminosity distance of 200 Mpc). The white con¬ 
tours show the projected gas pressure. 


feasible. In the following, we will explore the simulations we 
have presented here to extract as much as possible about cos¬ 
mic filaments in radio emission. We notice, however, that weak 
merger shocks leading to relics are more subject to resolution 
effects in grid methods than strong accretion shocks, and that 
a resolution of ~ 200 kpc is already sufficient to capture the 


bulk of energy dissipation in large-scale shocks (e.g. lRvu et al 


20031: 1 Vazza et al.ll2009il201ll) . All runs presented in this work 


have a spatial resolution better than this, and the estimated ra¬ 
dio power of our filaments is therefore a more robust statistics. 

Second, to better compare with the observed distribution 
of radio relics, we compute the radio distribution function for 
each magnetic field model in our intermediate (100 Mpc) 3 
volume at z = 0.05, and compare it with the relic distribu¬ 
tion function based on the NVSS and derived by iNuza et aD 
(120121) . This is so far the best guess for the abundance of radio 
relics in the local Universe, obtained by extrapolating the dis¬ 
tribution function of radio relics in simulations, with adjustable 
parameters calibrated using all radio relics in the NVSS and 
with X-ray scaling relations. Now we also introduce the ob¬ 
servational effect discussed in Sec 12.3.31 and mimic the sen¬ 
sitivity of the NVSS survey.The result is shown in Fig. [6) As 
expected, the baseline MHD model cannot reproduce the dis¬ 
tribution of bright radio relics, owing to its very low magnetic 
fields in clusters. The other models instead show a reasonable 
agreement with the statistics derived from NVSS, with sig¬ 
nificant gaps only at high radio power (> 10 26 W/Hz) due 
to the lack of massive objects, and and very low radio power 
(Fradio < 10 23 W/Hz) due to the combined effect of lim¬ 
ited spatial resolution to properly resolve the internal struc¬ 
ture of the smallest clusters in the box, and of the loss of low- 
surface brightness objects when the observational cuts are in¬ 
troduced. The differences between magnetic field models with 
post-processing amplification are very small because their dif¬ 
ferences only come into effect beyond the virial radii of clus¬ 
ters. 


Even if the DSA model for the weak shocks responsible for 
relics is still under debate, it is reassuring that our models do 
not overestimate the counts of radio relics in the NVSS. We 
assume that the extrapolation of these models to the stronger 
accretion shocks around filaments and the typical condition of 
the WHIM is valid. 


3.2. Which level of magnetisation can be detected? 

We start by focusing on the minimum magnetic field strength 
that can enable a detection of the cosmic web at any redshift. 
Fig. □ gives the total radio emission found with the HA 
and LA model for our most resolved run (50 3 Mpc 3 ) with 
a resolution of 20.8 kpc), with additional contours of gas 
pressure to trace the distribution of baryons. In both models 
the strongest emission patterns trace merger shocks inside 
clusters, producing structures similar to giant radio relics. The 
flux densities from these relics reaches ~ 0.01 — 0.1 Jy/arcsec 2 
in the HA model, ^ 10 -5 — 10 -4 Jy/arcsec 2 in the LA 
model. At the accretion shocks around clusters and filaments 
the differences are similar: the flux densities typically reach 
^ 10 -8 Jy/arcsec 2 in the HA model and ^ 10“ 11 Jy/arcsec 2 
in the LA model. No significant differences are found at the 
scale of filaments when the re-acceleration of electrons at 
shocks is neglected (not shown), which means that filaments 
only produce radio emission where M > 10. 
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Next we compute how much of the above emission patterns 
can be detected by real radio telescopes, and we start with one 
of the most promising cases explored here, the Large Survey 
with LOFAR-HBA (first configuration in Tab. 2, assuming a 
noise per beam of 0.25 mJy and a resolution beam of 25”) 
at 120 MHz (Fig. 0. The region covers a projected area of 
(50 Mpc ) 2 with ~ (205°) 2 at z = 0.05 and ~ (5.3 °) 2 at 
z = 0.5, while for comparison the field of view of a single 
LOFAR-HBA pointing is (100°) 2 . The comparison between 
the high and low amplification model shows that most bright 
regions inside of halos (considering projected spheres of ra¬ 
dius 1.5-Rvir at all epochs) are identified in both cases (blue 
regions in the figure), the detectable regions outside of ha¬ 
los strongly depend on the amplification model. In the HA 
model the LOFAR-HBA survey should detect several bright 
emission spots within a volume z < 0.5, in most of cases 
related to outer accretion regions surrounding galaxy clusters. 
The best chances of detecting ~ Mpc-sized shocked regions 
of the WHIM in filaments are for z ^ 0.1, which gives a 
good trade-off between having a good coverage of the neces¬ 
sary short baselines and not being too affected by the cosmo¬ 
logical dimming in surface brightness. In the LA model only 
sparse and more compact knots of bright enough emission are 
detectable, with a poor reconstruction of the underlying WHIM 
distribution. At 2 = 0.5 (and also up to z = 1, limited to a very 
few bright spots) only a magnetic field as high as in our HA 
model will allow some detection of gas not associated with the 
virial region of clusters. In the LA case, basically no detection 
will be possible. 


3.3. Performances of radio surveys 


Next we investigate in depth the chances of detecting the cos¬ 
mic web by existing and future surveys. First, we assess the 
role played by spatial resolution by comparing the maps of four 
low-frequency surveys, which are the most promising for a de¬ 
tection of the cosmic web in a high-amplification scenario (see 
also discussion below). Fig. [9] shows a 5° x 7° region for the 
highest resolution volume of (50 Mpc) 3 , observed by LOFAR- 
HBA at two resolutions, SKA1-LOW and MWA. The spatial 
resolutions here span from 10" (SKA1-LOW) to 120" (MWA). 
Each survey should be able to detect a few bright emission re¬ 
gions from the WHIM. However, only having a large enough 
resolution it is possible to guess the morphology and the spa¬ 
tial orientation of these features, which are otherwise difficult 
to be distinguished from the possible signal of relic-like emis¬ 
sion from distant clusters not yet detected in the X-rays (e.g. 
Macario et al.ll201 J) . Furthermore, achieving the highest pos¬ 


sible resolution is important for removing point sources and 
thereby for reducing the confusion noise, which is expected to 
limit the sensitivity of SKA1-LOW (e.g. Prandoni & Seymour 

12014 ) . 


Second, we study the performances of all surveys in Tab. 2 
using a more realistic setup, where we integrated a long rect¬ 
angular volume covering 14° x 14° in the sky and sampling all 
cosmic volume from z = 0.02 to z = 0.5 (corresponding to a 
comoving radial distance of « 1.892 Gpc). To do so we stacked 


many simulated volumes along the line of sight, starting with 
a few replicas of the (100 Mpc ) 3 box to give a better angular 
resolution at low redshift, and then by stacking several replicas 
of the (200 Mpc ) 3 volume to cover all projected volume up to 
z = 0.5. We first computed the radio emission in the comoving 
frame of each box at the appropriate redshift, and then “moved” 
each map away from the observer, by applying the cosmolog¬ 
ical corrections for the surface brightness and the luminosity 
distance of each redshift, and decrease the pixel size of each 
box as a function of distance, using a cubic interpolation on 
the input map. The presence of possible artefacts due to from 
the periodicity of structures along the line of sight is minimised 
by shifting each box with random offset in both directions. We 
notice that the first simulated volume along the line of sight 
is large enough that the change of redshift across the volume 
cannot be neglected. Therefore, we further subdivided it into 
4 slabs (starting at z & 0.023, 0.029, 0.035 and 0.044), and 
accordingly computed the cosmological corrections for each of 
them independently. 

Fig. [T0| shows effect of adding more volume along a pro¬ 
jected field-of-view, where we show the proxy of X-ray emis¬ 
sion (oc n 2 T 1 / 2 ) and the radio emission for a subset of the full 
map created above, assuming the HA model and for a Large 
Survey with LOFAR-LBA. Despite the evident appearance of 
many additional clusters in the field, when our integration re¬ 
gion goes from z = 0.03 to 2 : = 0.5 the detectable radio emis¬ 
sion does not change much, due to the strong cosmological 
dimming for 2 < 0.1. A few brighter spots, correlated with 
higher redshift merging clusters appear in the full-volume im¬ 
age, yet most of the cosmic web stays below the detectability 
threshold for LOFAR-LBA (large survey with resolution 25"). 

To give a more quantitative view of what fraction of the 
cosmic web can actually be detected by different surveys and 
for different magnetic field models, we give in Fig.[TT]the result 
for the full integrated volume. We show the total emission from 
all shocked cells, or from the WHIM only after excision of the 
projected virial volume of each identified cluster or group. In 
the HA model, the power in the WHIM is ^ 1% of the ra¬ 
dio power emitted by clusters, while this is ~ 0.1% in the LA 
model. Most of the emission in the projected volume comes 
from clusters and all surveys appear to be able to detect most 
of it (within a factor ~ 2 ), with the least detection coming from 
the highest frequency surveys. On average, we measure a dif¬ 
fuse radio flux density of r\j 0.1 Jy/deg 2 (100 MHz tv) from all 
shocks in our simulated volume (very similar numbers for both 
magnetic field model), while the contribution from the WHIM 
is ~ 5 • 10 -3 — 10 -2 Jy/deg 2 (100 MRz/v) in the HA case and 
- 10 - 3 - 5 • 10 “ 3 Jy/deg 2 (100 MUz/v) in the LA model. 

The differences in the performance between surveys be¬ 
comes more apparent if we look at the residual WHIM emis¬ 
sion only. The most satisfactory reconstruction of the WHIM is 
obtained with LOFAR (HBA Large 1 and LB A Large surveys), 
SUMSS, WENSS, SKAl-LOW, SKA1-MID and MWA-SUR, 
with a fraction of recovered flux which is 1/50 — 1/100 of the 
total WHIM emission. A bigger gap between the total and the 
detectable emission from the WHIM is measured for the LA 
case. 
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Fig. 8. Mock observation of the (50 Mpc) 3 volume using the configuration of LOFAR-HBA (Large Survey 1, see Tab. 12.3.21) for 
increasing redshift (z = 0.05, z = 0.1 and z = 0.5) for the HA (right) and LA (left) magnetic field models. The blue and red 
regions show the emission that is detectable for each model/redshift, divided into emission coming from < 1.5r v i r (blue) of each 
identified halo, or generated outside and therefore likely related to filaments (red). The additional contours show the projected 
gas pressure. The units of the axis are given in cell units (Ax = 20.8 kpc). 
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Fig. 9. Role of resolution in different radio surveys at low frequency. The colours show the detectable emission by different 
surveys at low-frequency for the HA model (units of log 10 [Jy/arcsec 2 ]), considering a 5° x 7° field of view. The emission in 
each map has been convolved with the corresponding resolution beam. The white contours show the projected gas pressure. 


3.4. Filaments connecting clusters 


Even before the completion of future deep radio surveys, it 
will be possible to detect the radio emitting WHIM with long 
exposures targeting galaxy clusters at low frequency. Close 
pairs of interacting clusters are particularly interesting for thi s 


(e.g. lPlanck Collaboration et~aD 201 31:iFarns worth et al.ll2Q13l) . 

Although these objects are much smaller than the ^ 10 — 
20 Mpc filaments that can be detected in our simulations, 
they can give us a first hint of the active mixing between the 
ICM and the WHIM. First DM-only simulations showed that 
most of simulated clusters are found to be surrounded by fila¬ 
ments with typical lengths < 20 Mpc/h, and with a weak mass- 
dependence: clusters with M ~ 10 14 M© are on average sur¬ 
rounded by ^ 2 — 3 filaments while this increases to ^ 4 — 5 
forM ~ 10 15 M© objects (e.g. lColberg et al.l200l) . Our recent 
analysis of the gas properties of the filaments in large cosmo¬ 
logical simulations showed that up to ~ 10 3 filaments can be 
identified in (100 Mpc) 3 at z = 0, with a decreasing power- 


law distribution of lengths from one Mpc up to several tens of 
Mpc (Gheller et al., submitted). Deep imaging of the region 
between close pairs of galaxy clusters is potentially an efficient 
way to search for cosmic filaments, even if in such close config¬ 
urations it might be difficult to disentangle the emission from 
cluster accretion shocks and filaments shocks. 

In a first exploratory study, we have extracted the position of 
halos with > 5 • 1 O 13 M 0 in our most resolved box and have 
extracted various quantities along the line joining the centres 
of close clusters. A more complete survey of the cluster-to- 
filament connections (involving the cross-correlation of clus¬ 
ter and filament catalogs) will be the subject of future work. 
Here we show a couple of representative cases where some 
emission from the WHIM can be detected. Fig. [12] shows the 
projected maps of magnetic fields (with isodensity contours) 
and of radio emission at 120 MHz (observed with FOFAR- 
HBA Farge Survey 1) for two pairs of interacting clusters in 
the (50 Mpc ) 3 box. The distances between the centres of the 
gas halos are ^ 4.5 and 6 Mpc, and the connecting bridges are 
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Fig. 10. Emission from our simulated light cone at z = 0.03, z = 0.05 and 2 = 0.5, obtained by stacking replicas of our 
100 Mpc) 3 and 200 Mpc) 3 simulated volumes. The top panels show the total X-ray emission proxy (oc n 2 T x / 2 ) and the lower 
panels show the radio emission for LOFAR-LBA (units of log 10 [Jy/arcsec 2 ]) for a projected area of 3.5° x 3.5°, assuming the 
high-amplification model (HA). The white contours of the X-ray emission proxy are drawn in the lower panels to allow a better 
comparison of structures. 


~ 1 Mpc thick. The magnetic field along the two filaments is 
~ 0.05 — 0.1/iG in the HA model and ~ 0.001/iG in the LA 
model (except that in a few dense substructures within the fila¬ 
ment), as shown by Figs. fl3lfl4l The interaction region of both 
pairs host significant radio emission, associated with several 
shocks within the filaments and around the gas substructures 
internal to it. In both cases for a high level of magnetisation 
a few bright emission regions ^ 100 — 200 kpc long can be 
detected, while most of the weaker shocks close to the spine 
of the bridges will stay below the sensitivity of LOFAR. In a 
low-magnetisation scenario basically no detection from these 
bridges will be possible, while the relic-like emission inside 
the main halos remains similar. The presence of radio galax¬ 
ies can add radio emitting plasma on these bridges, specially if 
such close pairs are found in a rich cosmic environment. 


3.5. Alternative scenarios 

We plan to investigate the role played by galaxies and ac¬ 
tive galactic nuclei in the seeding of large-scale magnetic 


fields with direct more complex numerical simulations in fu¬ 
ture work. Here we briefly discuss a few possible variations of 
the magnetic field models studied above, which we can easily 
test with existing runs. A possibility suggested in the litera¬ 
ture is that cosmic rays accelerated by strong shocks can in¬ 
duce turbulence in the shocks up-stream and trigger substantial 
magn etic field amplification (Drury & Downes 


2013). We tested this scenario (labelled HSA) by computing 


the kinetic energy flux which should be dissipated into CRs 
based on our shock finder and on the acceleration efficiency by 
Kang & Ryu ( 2013 ). This is converted into a magnetic field for 


each cell, by integrating for the cell size and the time step of 
the simulation: F> eq; cR = \/ 0 . 1 i?cR 87 r, i.e. we assume that a 
10% of the energy of accelerated CR goes into magnetic field 
amplification. The field obtained in this way is added to the 
field in the previous HA model. Fig. Q3] shows the radio emis¬ 
sion obtained in this scenario: the observed morphologies and 
power are very similar to the HA case, with the exception of 
the accretion shocks, where the magnetic field can be ampli¬ 
fied by a factor ~ 10 and the emission can be increased up to 
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B-model HA 



B-model LA 



Fig. 11. Performance of different radio surveys observing a 14° x 14° region. We simulated all radio emitting volume from 
z = 0.02 to 2 = 0.5 and present for both models the radio emission from the ICM+WHIM (blue) or from the WHIM only (red) 
distinguishing the total emission (big squares) and the one detected by each survey configuration (small diamonds). For a better 
comparison of surveys at different frequencies, all emission have been rescaled by (110 MHz/i/[MHz]). 


^ 10 — 100 at accretion shocks around clusters and filaments. 
However, this extra amplification is often found at densities too 
low to contribute significantly to the observable radio emission. 
Fig. Q2I shows the outcome of the HSA model on the emitted 
and detected radio emission for the same long projected vol¬ 
ume analyzed in the previous sections. The total emission in 
the box, largely dominated by clusters, is the same when the 
HA and HSA model are compared, and also the total emission 
from the WHIM remains basically the same. In any case, for 
the moment this mechanism remains just a speculation, given 
the largely unknown acceleration efficiency of CR-protons by 
cosmological shocks, and the complex details related to the dif¬ 


fu sion coefficient of particles in the upstre am (see discussions 
in iDrurv & Downesll2Q12l:lBruggenll2013h . 


Secondly, we tested the case in which the radio emitting 
electrons are only those freshly accelerated by shocks and there 
is no contribution from re-accelerated electrons at low Mach 
number (Sec. 12.3.21) . Fig. Q2I shows that neglecting shock re¬ 
acceleration decreases the total emission by a factor ~ 2, 
by reducing the relic emission in clusters (as already seen in 
Sec. 13.1b . Also the WHIM emission is reduced by a similar 
factor. 


In summary, the presence of additional amplification by 
cosmic-ray driven turbulence or the presence of pre-existing 
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Fig. 12. Maps of intracluster filaments in the 50 Mpc) 3 volume 
at z = 0.05. The left panels show the projected mass-weighted 
magnetic field (HA model, units of log 10 [/iG]) and the white 
contours give the projected gas density. The right panels show 
the intrinsic radio emission in the same model ([Jy/arcsec 2 ]). 


electrons should affect our estimates only by a very small fac¬ 
tor, < 2 in the total and detectable emission. 


Finally, we notice that the total emitted power from the 
WHIM essentially scales linearly with the saturated value of 
£(M, T) for strong shocks (Eq.Q]). In this work this is £ ~ 0.7— 
0.9 • 10 -4 , of the same order of what usually assumed to ex¬ 
plain the observed spectra of supernova remnants dEdmon et al. 
201 ltlFerrand et all2014h . Given that the predicted radio spec¬ 
tra of most shocks in the WHIM is flat P v oc z/ -1 , Eq. Q] pre¬ 
dicts that for magnetic field lower than 3.2 fiG and at low red- 
shift, different combinations of £ • B 2 should give the same 
emission level. Considering the level of predicted emission 
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Fig. 13. Profiles along the white line in the first panel of 
Fig.UJl showing the trend of mean gas density along the line of 
sight, magnetic field (each color gives a different model, with 
blue=LA and green=HA), intrinsic radio emission (same color 
coding) and detectable radio emission assuming the sensitivity 
of LOFAR HBA Large Survey at 120 MHz (resolution 25”). 


from the WHIM in our runs (Sec. 13.3b . we can therefore re¬ 
formulate our model prediction in the following way: 

-FVhim(z') ~5x 10 -3 Jy/deg 2 -—-- pjp • ppj , (2) 

where cb is the magnetisation level of the WHIM we as¬ 
sumed in the HA model, i.e. 1% of the thermal gas energy. 
This can also be roughly translated into a prescription for the 
mean magnetic field of the WHIM, assuming a mean density 

of n/n CY ~ 10: 

PwhimM ~ 5xi°- 3 . - 1 ^.( 3 ) 

Thus, the level of radio emission we predict for the WHIM can 
be rescaled for the simple cases of a different constant magneti¬ 
sation level of the WHIM (e.g. assuming a different cosmolog¬ 
ical seed magnetic field) or of a revised acceleration efficiency 
for electrons (e.g. revised acceleration efficiency from super¬ 
novae). 

4. Conclusions 

In this work we studied the detectability of the cosmic 
web using existing and future radio surveys. We simulated 
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Fig. 14. Same as Fig.[l3l but for the second cluster association 
in the Figure. 
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Fig. 15. Similar to Fig. [71 but using the HSA model (high am¬ 
plification and CR-driven dynamo). 


SKA 1 — LOW, different models 



Fig. 16. Comparison of the emitted and detected radio emission 
for the same 14° x 14° region of Fig. Qj] using SKA1-LOW. 
Each column shows the outcome of different models for the 
magnetic fields (HA, LA, HSA) and also without the contribu¬ 
tion from re-accelerated electrons. 


the acceleration of elec trons by cosmological shock waves 
( Hoeft & Bruggenl 12007 ) and applied it to large cosmolog¬ 
ical volumes produced using the MHD version of ENZO 
d Vazza et al.ll2014ah . In post-processing we modified the mag¬ 
netic field as a function of density such as to reproduce the 
magnetic filed levels observed in galaxy clusters. We increased 
the magnetic field according to a high-amplification (HA) or 
low-amplification (LA) models, which both pushes the mag¬ 
netic field in clusters close to the observed values (a few ^ fiG) 
but differ in filaments: in the HA model the magnetic field in 
filament is set to 1 % of the thermal gas energy of cells, while 
in the LA it remains the one from the MHD simulation, i.e. 
< nG. We produced mock radio observation with the specific 
parameters of real surveys (including the lack of short baseline, 
of the spatial resolution and sensitivity of each instrument) and 
quantified the chance of detecting the WHIM in large projected 
volumes (up to z = 0.5). 

Our main results can be summarised as follows: 


- Only the high-amplification scenario allows a detection of 
the tip of the iceberg of the WHIM radio emission, lim¬ 
ited to the brightest emission spots in filaments and accre¬ 
tion shocks. The high-magnetisation might come from un¬ 
resolved small-scale dynamo (or fast growing instabilities 
in a high plasma /3) and/or from seeding mechanisms not 
yet included in our MHD runs (e.g. outflows from galaxies, 
AGN), which push the magnetic field up to ~ 1% of the 
thermal energy of the WHIM as in clusters (Sec. i m 

- Low-frequency (< 300 MHz) observations are best suited 
to detect the large-scale diffuse emission produced by cos¬ 
mological shocks, at the scale of cluster accretion shocks 
or in filaments. These shocks are usually several ~ Mpc 
wide, are characterised by a flat emission spectrum (P v oc 

1_ ^ 1 • 5 ) and emit at the level of ~ /iJy/arcsec 2 at low 
redshift, in a high amplification model (Sec. 13.3b . 

- Most of the detectable radio emission is caused by struc¬ 
tures at z < 0.1, because of cosmological dimming. 
Especially at high frequencies, the detection of too nearby 
structures of the cosmic web (z < 0.02) is made difficult 
by the lack of short baselines. 
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- Of all investigated surveys, the highest fraction of detected 
flux from the WHIM phase is obtained by LOFAR (HBA- 
Large survey and LBA-Large), SKA1-LOW and MWA 
Broadband survey. In all these cases 10“ 2 of the WHIM 
emission can be detected, mostly associated with shocks in¬ 
ternal to filaments, both driven by the accretion of smooth 
gas in rich environment or driven by supersonic motions 
internal to filaments. 

- The highest resolution accessible to LOFAR (< 25”) and 
SKA1-LOW (~ 10”) should enable a much better identi¬ 
fication of the morphology of the underlying filamentary 
environment, compared to MWA, and to exclude the con¬ 
tamination from cluster shocks and active galaxies. High 
resolution will be also important to reduce the confusion 
noise (e.g. in the SKA1-LOW). 

- In a low-amplification scenario the detection chance of 
the above surveys will be reduced by ^ 10, and only a 
very limited number of bright spots (mostly associated with 
dense substructures within filaments) will be detectable, of¬ 
fering a poor tracing of the WHIM distribution. 

- We estimate a background of unresolved radio emis¬ 
sion from the WHIM in the range ^ 10 -3 — 

10 -2 Jy/deg 2 (100MHz/V), depending on the amplifica¬ 
tion model (Eq.03]). 

- Clusters interacting on short distances 4 — 6 Mpc) can be 

connected by magnetised bridges, hosting internal shocks 
that can be observed if the magnetisation level is ~ 1% of 
the thermal gas (Sec. 13.41) . 

- Our modeling does not yet cover the impact of magnetised 
outflows from galaxies and active galactic nuclei in a proper 
way. We investigated a few alternative scenarios for the am¬ 
plification of magnetic field and for the acceleration of par¬ 
ticles at strong shocks, finding only minor differences with 
respect to the main results of our baseline model (Sec. I331> . 

- Radio emission from secondary electrons injected by 
hadronic collisions is negligible in the WHIM (Appendix). 


An important caveat in our analysis is that we assume that 
DSA can operate at all density and magnetisation levels in our 
simulated volume, which may not be true for very low levels 
of magnetic fields. Recent results from Particle In Cell simu¬ 
lations (PIC) pointed out that more important for the onset of 
DSA at non-relativistic shock is the topology of the upstream 
magnetic field, rather than its initial value. In particular, if the 
upstream magnetic field is aligned with the shock normal (i.e. 
parallel shocks), the rapid development of the Weibel insta¬ 
bility might amplify the magnetic field and enable^ the _for- 
mation of a shock precursor leading to DSA (ICapriolil 12012 : 
ICaprioli & Spitkovskvll20l4lGuo et al .1120141) . For this reason, 
the presence of shock-accelerated electrons in the WHIM is a 
realistic possibility. 

An important aspect of the above scenario is that both a de¬ 
tection or a non-detection of the cosmic web in radio within 
the next decade will be helpful to put constraints on the mag¬ 
netic field in the WHIM, given that to a first approximation the 
radio emission from the WHIM should scale as P v oc £ B 2 
(Eq]2]-Eql3]). The lack of detection in radio from large ob¬ 
served volume will limit the magnetic energy in the WHIM 


to values much below ~ 1 % of the thermal gas energy in 
these environment. In this case, the effect of cosmic magnetic 
fields on the propagation of UHECRs is p redicted to be negli¬ 
gible dSigl et al] 20031: iDolag et aPl2005l) . and so is the effect 
of magnetic fields in producing axion-like particle oscillations 
dHorns et aUj201 2). 


On the other hand, in the case of systematic detection of 
emission patches connected to the WHIM, the analysis of their 
morphology and large-scale distribution will be important to 
distinguish competing scenarios. For example the detection of 
coherent magnetic fields on >> Mpc scales will favour a dis¬ 
tributed amplification process as an obliquitous small-scale dy¬ 
namo. Our models show that if the growth of cosmic web is 
linked to the thermal gas energy of the WHIM, the cosmic web 
should be maximally emitting in radio at moderately low red- 
shift (z < 0.1). Any unambiguous detection of radio emission 
from much higher redshift will further discriminate between 
seeding models and favour a strong magnetisation activity of 
large-scale structures by high-2 star-forming galaxies. 


Beside continuum radio emission, other complementary 
techniques should be able to probe the existence of the cos¬ 
mic web, using differe nt effects: polarised r adio emission 
(e.g. Govoni et aPl2015h, Faraday Rotation (e. g. lAkahori et al.1 
2014 : lBonafede et al .1120151: iTavlor et al.ll2015l) . neutral hydro¬ 


gen emission (e.g. Popping et al. 2015|), absorption and emis- 
sion features of hig hly ionised elements i n the WHIM (e.g . 

and fast radio bursts (Zh eng et al] 12014 


Barcons et al. 201 


Dolag et al.1120141) . Considering the extreme difficulty of get¬ 


ting the expected weak signal from the cosmic web, the use of 
statistical an d stacking techniques will increase the chances of 
success (e.g. IVacca et‘a.1 J201.4lSt.il & Keltefeoi 5h . 
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Appendix A: The contribution from secondary 
electrons 

We briefly report here on the contribution from sec¬ 
ondary electrons injected by hadronic proton-proton collision 
(Blasi & Colafrancescol l 19991: Dolag & EnBlinl EoOOh to the to¬ 
tal radio emission. The injection rate of secondary electrons 
is computed from the pool of CR-protons simulated at run¬ 
time by our ENZO cosmological simulation using PPM and a 


two-fluid model for modeling cosmic rays, as in Vazza et al 


( 2014 b). For this we use the formalism by Dolag & EnBlin 
<I2000L assuming a constant energy spectrum of a = 1.1 
(. N(E ) oc E ~ a ) for cosmic ray protons. The final bud¬ 


get of cosmic rays in the simulation follows from the as¬ 
sumed acceleration efficiency at shocks, which we took from 
Kang & Rvul ( 2013 ) and renormalize downwards of a factor 10 


to be consistent with the upper limits from the Fermi satellite 
dAckermann et al.l 12014 ). Here we use the result from a (300 
Mpc) 3 volume simulated with 2048 3 cells and DM-particles. 
This run makes does not use MHD and therefore we have to 
assume the magnetic field entirely in post-processing. Similar 
to the main articles, we compute the gas energy in all cells and 
compute the magnetic field in a HA and FA model. 

Fig. IA.1I shows the outcome of these two acceleration 
mechanisms for the HA model: the contribution from sec¬ 
ondary electrons to the diffuse cosmic web is everywhere neg¬ 
ligible, i.e. only a few percent of the primary contribution. The 
emission gets stronger in the centre of the clusters owing to 
the increase with gas densit y, an d only forms rather low-power 
radio-halos (e.g. Dolag & Enfilinl I2OO0I) . Outside the central 
cluster regions, however, this signal usually gets much smaller 
than the primary emission from strong accretion shocks. The 
distribution function of pixels in the maps (Fig. IA.2I) shows 
that by far the primary contribution dominates the high-end of 
the emission tail in both magnetic field model. However, the 
secondary contribution becomes the dominant one in the low¬ 
brightness end of the emission distribution in a FA scenario, 
due to drop of primary emission from the WHIM outside of 
halos and to the fact that the secondary emission is nearly un¬ 
altered in the FA scenario since it mostly comes from well 
within halos. The confusion caused by the secondary emis¬ 
sion from galaxy clusters should become more important at 
high redshift, where massive bright clusters can dominate the 
low power background emission compared to the filaments, in 
the FA scenario. Overall, we conclude that regardless of the 
assumed amplification model, the brightest emission from the 
cosmic web should be originated by primary electrons acceler¬ 
ated by cosmic shocks, and that the signal from secondary par¬ 
ticles can be relevant only at very low values of surface bright¬ 
ness and at high redshift, in a low magnetisation scenario of the 


WHIM. 
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Fig.A.2. Total radio emission from primary electrons injected by shocks (right) and from secondary electrons injected by 
hadronic collisions (left) in a subvolume of our 2048 3 simulation with CR-physics dVazza et al.ll2014bl) . The color bar gives 
the emission in in units of [logioW /Hz]. 



F[Jy/arcsec 2 ] 


Fig. A.l. Distribution of radio emission from the HA and the 
LA model in a cosmological run with CR-protons and includ¬ 
ing the effect of “secondary” radio emission from electrons in¬ 
jected by hadronic collision. 











